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Abstract 

In systems undergoing second order phase transitions, the temperature integral of the 
specific heat over temperature from zero to the critical temperature is the same in both 
the normal and ordered phases. This entropy balance relates the critical temperature 
to the distribution of degrees of freedom in the normal and ordered states. Quantum 
criticality and fractionalization can imply an increased number of low energy degrees of 
freedom in both the normal and ordered states. We explore the role of entropy balance 
in holographic models of superconductivity focussing on the interplay between quantum 
criticality and superconductivity. We consider models with and without a ground state 
entropy density in the normal phase; the latter models are a new class of holographic 
superconductors. We explain how a normal phase entropy density manifests itself in 
the stable superconducting phase. 



1 Introduction 



Perhaps the central question concerning superconductivity is the most obvious one: when 
can the critical superconducting temperature become high? What circumstances favor 
quantum mechanical pairing instabilities against the stabilizing effects of thermal fluctu- 
ations? The theorist might like to identify aspects of these big questions that can be 
abstracted as much as possible from the specifics of concrete materials. 

As we will shortly recall, systems undergoing second order phase transitions obey a 
type of entropy balance: the integral of the specific heat over temperature, from zero 
temperature to the critical temperature, has to be the same when evaluated in the (stable) 
superconducting state and in the (unstable) normal state. This gives a formula that relates 
the critical temperature to the difference in the distribution of degrees of freedom between 
the normal and superconducting states. Read in the most naive way, it suggests that having 
high critical temperatures is correlated with the normal state having many more low energy 
degrees of freedom than the superconducting state. For conventional superconductors the 
entropy balance is well understood, as we will briefly review: the superconducting state 
is gapped, the number of 'surplus' low energy degrees of freedom in the normal state are 
proportional to this gap, and consequently T c is also proportional to the gap. 

Unconventional superconductors have unconventional normal states. There is broad 
evidence for the notion that the normal states are unconventional because the low energy 
fermionic degrees of freedom are strongly interacting with collective critical bosonic exci- 
tations. See e.g. [1] for a suggestive overview of several different compounds. In terms of 
the entropy balance, additional low energy degrees of freedom in the normal state bodes 
well for obtaining a higher critical temperature. It is important to have concrete, controlled 
models where the entropy balance can be explicitly worked out in such cases. For instance, 
as we shall see, it can occur that the critical degrees of freedom are not all gapped out in 
the superconducting state. In this paper we shall see how the entropy balance pans out in 
holographic models of superconductivity. As we shall also briefly review, these are theo- 
retical models for the emergence of superconductivity from fractionalized phases of matter, 
wherein the charge carriers are strongly coupled to deconfined gauge fields. 

An important impetus for our study were the recent careful measurements of the tem- 
perature dependence of the specific heat in Sr3Ru207 [2]. The evidence suggests that, in a 
suitably tuned magnetic field, this material is a two (space) dimensional quantum critical 
metal that undergoes nematic ordering at low temperatures, as reviewed in [3]. The general 
considerations we have described above apply equally to nematic as well as superconducting 
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instabilities. Among the several dramatic discoveries in [2] is that, as the temperature is 
lowered through the critical nematic temperature, the specific heat over temperature contin- 
ues to increase in the nematic phase as the system is cooled. It seems natural to think that 
this increase is made possible by the persistence of quantum critical degrees of freedom in 
the ordered state. In the models studied in this paper, which are presumably not related to 
Sr3Ru207, we will precisely explore the effects of quantum criticality on the redistribution 
of degrees of freedom as the system undergoes an ordering transition. 

A second motivation for this work is that, as we will recall below, many of the simplest 
models of holographic superconductivity generically have the property that the normal 
state of the system is locally quantum critical and exhibits a zero temperature ground state 
entropy density. This state of affairs is usually considered to be 'cured' by the onset of super- 
conductivity in the sense that the superconducting state at zero temperature has vanishing 
entropy density. Nonetheless, we will see that entropy balance implies that the normal and 
superconducting states are not independent; the zero temperature entropy density of the 
normal state can have important consequences for the thermodynamics of the supercon- 
ducting state. For comparison, we will study holographic models of superconductivity with 
and without a ground state entropy in the normal state. 

Finally, we have also been inspired by recent applications of the notion of "quantum 
order-by-disorder" to finite density systems, namely, Fermi liquids in the vicinity of ferro- 
magnetic instabilities [4, 5, 6]. The idea is that the low energy degrees of freedom associated 
with disordering the original ferromagnetic phase can, as the quantum critical point is ap- 
proached, favor new orderings. Technically this is described in [4, 5, 6] using second order 
perturbation theory about a mean field treatment. One would like to have a general model- 
independent formulation, and furthermore not tied to a perturbative understanding of the 
difference between classical and quantum effects, of what quantum order-by-disorder means. 
The notion of entropy balance that we will shortly review in detail has the flavor of a gen- 
eralized notion of quantum order-by-disorder, as it can capture general features relating to 
the interplay of quantum criticality and ordering. 

2 Second order phase transitions and entropy balance 

At a second order phase transition, the free energy and entropy densities of the supercon- 
ducting and normal states are both equal at the critical temperature 

/jv(T c ) = f s (T c ) , s N (T c ) = s s (T c ) . (1) 
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The specific heat, in contrast, will jump discontinuously at the critical temperature and, 
just below T c , will be larger in the superconducting state 



This inequality expresses the fact that the superconducting state is preferred for T < T c . 
The above statements are strictly true if the transition has a mean field form. This will be 
the case for all of the models considered in this paper. 

From the equality of the entropies at the critical temperature, and from the definition of 
the specific heat as c/T = ds/dT, one obtains the following 'sum rule' relating the integrals 
of the specific heat in the superconducting and normal phases 



On the left hand side we have allowed for a finite entropy density at zero temperature in 
the normal state, as this will be present for some of the models we consider. The relation 
(3) imposes an 'entropy balance' constraint that keeps track of the degrees of freedom, even 
though the excitations of the system can be drastically re-organized due to the onset of 
superconductivity. 

The entropy balance equation (3) suggests the following observation. Suppose we wish 
to tweak the system to increase the critical temperature. How might we do it? The in- 
equality (2) tells us that the contribution from near the upper limit of the integrals, just 
below T c , is larger in the superconducting state. Therefore, if we increase T c , the right 
hand side of equation (3) becomes larger. To maintain the entropy balance it follows that 
at temperatures well below the critical temperature we might either (a) increase the low 
temperature specific heat of the normal state or (b) decrease the low temperature specific 
heat of the superconducting state. We could of course do both. A, somewhat naive, physi- 
cal conclusion would seem to be that to increase the critical temperature we should try to 
increase the low energy degrees of freedom of the normal state relative to the low energy 
degrees of freedom of the superconducting state. An alternative way to increase the critical 
temperature is to decrease the jump in the heat capacity (2) at T c . This paper will be an 
investigation of such correlations in holographic models of superconductivity. 



To contextualize, we will very briefly review how the entropy balance works out in the BCS 
theory of superconductivity. We will only consider a mean field treatment here. The key 



c N (T c ) < c s (T c ) . 



(2) 





3 BCS theory 
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formula is of course the famous gap equation 

1 r D de , E 



, „ tanh — , E = e + Ai . (4) 
gN(0) J E 2T' T v 1 

Here g is the electron-phonon coupling and iV(0) the density of states at the Fermi surface. 
This equation immediately gives the zero temperature gap and the critical temperature in 
terms of the gap 

A = 2uj D e~ 1 /9 N ^ , T c « 0.57A . (5) 
The specific heat of the superconducting and normal phases are (e.g. [7]) 

c s = N(0) tfesech 2 — • — --— , CJV = — iV(0)T. (6) 



.,,„ 2T V r2 T dT J ' 3 

From this formula we obtain the low temperature specific heat in the superconducting phase 



/8^A^(0) (A /T) 5 / 2 e" A °/ T + • • • . (7) 

Inspecting equations (5) and (7) we see that the naive intuition from above is realized: 
Increasing Ao increases T c and simultaneously decreases the degrees of freedom in the su- 
perconducting state at low temperatures. Plotting the specific heat of the two phases in 
figure 1, for small gN(0), we can see exactly how the entropy balance works out. 
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Figure 1: Entropy balance in BCS theory. Specific heat over temperature for the normal 
and superconducting states below T c . The area between the two curves integrates to zero. 

The main feature in the BCS specific heat plot is the gap in the superconducting phase at 
low temperatures. In the holographic models that we will be studying in the following, the 
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normal state is quantum critical and not all of the quantum critical degrees of freedom will 
be gapped in the superconducting state. We wish to use holographic models to examine the 
interplay between criticality in the normal and ordered states and the critical temperature. 

4 Holographic superconductors 

4.1 The condensed matter physics of holographic superconductors 

Holographic superconductors are a gravitationally dual description of the onset of super- 
conductivity from certain fractionalized phases of matter. By fractionalized phase we mean 
that, somewhat heuristically, the electric charge may be thought of as being carried by 
degrees of freedom that are charged under an emergent large N gauge symmetry. These 
could be, for instance, N x N matrix-valued fermion fields \f or boson fields <£. This 
gauge symmetry is in addition to the global U(l) 'electric' charge. The superconducting 
instability occurs when it becomes favorable for the charge to be carried by short gauge- 
invariant bosonic operators, such as Tr (^ f \E f ) or Tr ($$). We will need two different flavors 
of fermionic field in order to make the trace of the product nonvanishing. These are not 
Cooper pairs of the UV gauge invariant 'physical' electrons, but rather pairs of low en- 
ergy fractionalized excitations. The physical electrons themselves are given by operators 
of the schematic form Tr(^<l>). The corresponding more conventional Cooper pairs would 
then correspond to the 'double trace' operator Tr (\]/<J>) Tr Such pairing can also be 

described holographically [8]. 

The fractionalized phase is gravitationally described by a charged black hole horizon. 
The formation of the 'meson' Cooper pairs described above is captured by the expulsion 
of charge from the horizon into a charged Bose-Einstein condensate outside the black hole. 
In the zero temperature limit, often all of the charge is expelled, yet a (neutral) horizon 
remains. This indicates that electrically neutral critical excitations, such as deconfined 
gauge fields, have survived the onset of superconductivity. This phenomenon is reminis- 
cent of the Landau damping of transverse gauge fields by a charge density, in which the 
transverse modes are not screened. The presence of critical modes in both the normal and 
superconducting states at low temperatures will be a recurring theme in this paper. 

A more detailed discussion of the above ideas appears in [9], along with references to 
the original literature. 

The discussion above shows that two processes occur simultaneously at T c in holographic 
superconductors. The global U(l) symmetry is spontaneously broken by a condensate and 
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furthermore the fractionalized charge starts to 'confine' into meson-like degrees of freedom. 
This is a little different to the emergence of superconductivity from a fractionalized phase 
in e.g. [10], where the analogue of the gauge-charged boson <3? condenses. In that case the 
global U(l) is spontaneously broken and the emergent gauge symmetry is Higgsed. Current 
holographic superconductor models do not appear to involve Higgsing of the SU (N) gauge 
symmetry. Pairing analogous to Tr (\E n I') was not accessible in [10] because the fractionalized 
fermionic excitations were electrically neutral. Other models of phases with emergent gauge 
fields such as the 'doublon metal' do have pairing instabilities in which the condensate is 
gauge- neutral but electromagnetically charged [11]. Various theories with similar phases, 
involving the simultaneous onset of confinment in the charged sector and superconductivity, 
are described in [12]. These are very much analogous to holographic superconductivity. 

4.2 Models with a normal state entropy density at T = 

We will proceed to study two classes of holographic superconductors. The first, simplest 
class, will have the curious feature that the normal state is locally quantum critical and 
has a zero temperature entropy density. Perhaps unsurprisingly, we will see this fact, the 
zero temperature entropy density, strongly colors the way in which the entropy balance is 
realized in these theories. In the following subsection we will move on to models that display 
a more generic quantum criticality in their normal states. 

A simple class of models of holographic superconductors have the bulk Lagrangian 

£ = ^ (r + ^) - - ^ (IW - iM? + m 2 M 2 + v(\d>\)) . (8) 

We will work with a 3+1 dimensional bulk, and hence a 2+1 dimensional dual quantum 
field theory. We will very briefly review the properties of the normal state of such models, 
with = 0. A more pedagogical exposition with the same conventions as we are using 
here can be found in [9]. The normal state is holographically described by the planar 
Reissner-Nordstrom-AdS black hole, with metric and Maxwell potential 

ds 2 = ^ (-f{r)dt 2 + ^+dx 2 + dy 2 ^j , A = M (l - dt . (9) 

Here \i will be the chemical potential of the dual field theory. The metric function is 

We introduced the dimensionless ratio of the Newtonian and Maxwell couplings 



The horizon radius r + determines the temperature as 



T =*k{*-%)- (12) 

The entropy density and specific heat of the solution are 

_L 2 2vr _ rr ds _ L 2 1 8vr 2 rr + 

SN ~^r{' CN - 1 df-^^3-2nTr + - (U) 

The temperature derivative has been taken with [i fixed. 

This class of models has the property that at zero temperature, the near horizon geome- 
try of the extremal Reissner-Nordstrom-AdS black hole is AdSi x M? . This implies that the 
low energy physics develops an emergent local quantum criticality, with scaling symmetry 
{t — > At, x — > x} in which space does not scale. Associated with this symmetry is the widely 
discussed zero temperature entropy density 

%(0) = ^- (14) 

Such an entropy density is not a necessary consequence of local criticality [13] , but does 
arise robustly in holographic models described at zero temperature by a regular extremal 
event horizon. 

In these models, the Reissner-Nordstrom solution can become unstable towards con- 
densation of 4> below a critical temperature T c [14, 15, 16]. The critical temperature itself 
depends on the mass m 2 of the scalar field as well as ratio 7 [17]. Despite the strong inter- 
actions of the dual quantum field theory, in which operators certainly have large anomalous 
dimensions, the finite temperature phase transition is mean field in nature because the large 
N limit suppresses fluctuations [18]. To discuss the superconducting state, we need to make 
a choice for the scalar potential and mass. 

A convenient set of models to consider have [19] 

m 2 L 2 = -2, vm = ^\ A - (15) 

This choice of the mass corresponds to taking the operator O dual to <fi to have scaling 

dimension A = 2 in the conformally invariant UV fixed point theory. The choice of potential 

ensures that at zero temperature in the superconducting phase the scalar field does not run 

off to infinity in the IR geometry. 1 This will make it easier to determine the low temperature 

1 While the slightly simpler case of positive mass squared and no potential also admits superconducting 
Lifshitz solutions in the IR at zero temperature [19, 20], with a stabilized scalar, fluctuations about these 
backgrounds always include modes with complex scaling dimensions and hence the geometries are unstable 
[21, 22]. Presumably, the true IR geometry in these cases has a runaway scalar. 
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specific heat of the system. In fact, the near horizon geometry at zero temperature has a 
scaling invariance characterized by a dynamical critical exponent z 

l l2 dt 2 dr 2 dx 2 + dy 2 . dt , 

72 ds = 27 + # L_ 2~ H 2 ' ^ = 7«L— , 4> = <PL- (16) 

Here /it,, </>.t} are constants. The exponent and stabilized value of the scalar field satisfy 

7 2 0|(4 + z + z 2 ) + z(^l-4</)|-6) = 0, (17) 
67 2 0i + z(u0!(3 + z)-24(4 + z)-6) = 0. (18) 

Furthermore h\ = (z — l)/z ,gL = z/(-f 2 <fi 2 L ). Note that z > 1. 

The 'Lifshitz' [23] scaling of the near horizon metric (16) at zero temperature implies 
that at low temperatures 

c s (T) ~ T 2 /* . (19) 

This follows from dimensional analysis. Thus the larger z is, the more degrees of freedom 
are present in the superconducting phase at low temperatures. Because z is a directly 
physical quantity, we will parametrize the theories by {z, 7} rather than {u, 7}. Given z 
and 7 it is easy to obtain the corresponding u from (17) and (18). The low temperature 
specific heat (19) is qualitatively different from the gapped BCS formula (7). As discussed 
above, it indicates that the critical modes in the normal state are not all gapped by the 
superconducting condensate in these models. 

The full spacetime, at zero and finite temperature, takes the general form 

1-ds 2 = -f{r)dt 2 + g{r)dr 2 + — % ^ , A = 1 h(r)dt, <f> = 4>{r) . (20) 

In these coordinates r — > will be the AdS± conformal boundary. This ansatz leads to the 
following equations of motion 

r^- ^ + 9 - + 6 - r 5 (6 + 40 2 - n^ 4 ) = 0, (22) 
f f 9 r 

h " + £ + rfghcf? = , (23) 



r 2/ 2g 

r 2f 



+ Q^ + {2g + 1 29 -^4>-ug^ = 0. (24) 



The strategy to obtain the specific heat as a function of temperature is now as follows. We 
must first find a solution to the above equations such that near the boundary <p ~ (0)r 2 
is purely normalizable. Away from the boundary, the solution should end at a regular 
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black hole horizon at r = r + . Near the boundary we can read off the speed of light from 
/ ~ c 2 r -2 . This is important because we must rescale the time coordinate so that all 
theories are compared with the same speed of light. By varying data at the horizon, we will 
find a one parameter family of solutions, from which we can obtain the specific heat using 
the following expressions for the temperature and entropy density 

1/2 L 2 2vr 



1 

47TC„ 



dfdig- 1 ) 



dr dr 



• ^ " K 2 r 2 • (25) 

r=r + K r + 



The family of solutions will correspond to temperatures ranging from T = to T = T c . 
It is important to compute dimensionless quantities; in fact we should obtain c/ (/xT) as a 
function of Tj ji, where the chemical potential is read off from the near boundary behaviour 
of the Maxwell potential: h ~ c fi/~/. We are imagining varying the temperature with the 
chemical potential held fixed. 

As usual, the equations are solved numerically by integrating out from the horizon and 
shooting to impose the boundary condition (f) ~ (O)r 2 . In general multiple families of 
solutions will be found and it is important to take the stable one, with the highest T c . As 
noted in [19], the existence of a Lifshitz scaling solution (16) does not guarantee that this 
solution arises as the IR limit of an asymptotically AdS spacetime. The condition for this 
to occur is found by perturbing the solutions away from Lifshitz 

/ = ^( 1 + /i ra )' 9 = 9 -±{i + gir a ) , h = ^(l + h ir a ), = ^(1 + ^0, 

(26) 

and requiring that there are at least two values of a < 0, such that the perturbations are 
irrelevant, i.e. grow towards the UV, and can be tuned. Substituting the above expansion 
into the equations (21) we will get five allowed nonzero exponents a; three are always 
positive and the other two can be negative depending on the value of 7. See [19] for a plot. 
In this way, a portion of possible values of {z, 7} do not in fact arise. We also require that 
a be real, as complex a indicates an instability [21, 22]. 

We saw in (19) that the number of degrees of freedom at the lowest temperatures in 
the ordered phase is determined by z. However, the critical temperature T c in fact only 
depends on 7 and not on z at all. We see this by recalling that T c is determined by a 
linearized analysis, where the charged scalar field is small. In this limit the nonlinear terms 
in the potential V(|</>|) in (8) drop out, and so, at fixed mass m, the equations of motion for 
the scalar field only depend on 7. The critical temperatures as a function of 7 for a range 
of masses were found in [17], using slightly different conventions to the present paper. In 
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particular it was shown that T c j ' \i increased with increasing 7, as this may be thought of as 
the effective charge of the scalar field. 

The independence of T c and z therefore seems to contradict our intuition that, if the 
normal state is kept fixed, there should be an anti-correlation between the critical tem- 
perature and the number of low energy degrees of freedom in the superconducting state. 
However, this is not the whole story. In figure 2 below we plot c/T against T < T c for the 
relatively small value of 7 = 1, and for three values of z = 1.5, 2, 3. 2 The first important 
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Figure 2: Specific heat over temperature as a function of temperature for the normal (red) 
and superconducting (blue) phases below T c and with 7 = 1. The plots are for three different 
values of z, for a model with a nonzero ground state entropy density Sjv(0). The shaded 
region between the blue curve and the dashed yellow line, indicating cjv/T + sn(0)/T c , 
integrates to zero area. 

feature of these three plots is that they are indistinguishable. We knew that T c would be 

the same, as it only depends on 7, but the low temperature behavior should have depended 

strongly on z; from (19) we see that cg/T goes to zero for z = 1.5, to a constant for z = 2, 

and diverges for z = 3. 

In this paragraph and the following, we explain how the plots in figure 2 can be nicely 

understood from entropy balance, despite being dramatically different from the BCS curve 

of figure 1. We see that for this value of 7, T c //j, is small, of order 10 -2 . Correspondingly, 

2 The Lifshitz solution with 7 = 1 and z = 3 is not in fact the true IR in this case, as there are operators 
with complex scaling dimension, i.e. a in (26), indicating an instability. This is not crucial for the points 
we are making here, and will not be the case in any other solutions we consider. 
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we also see that the superconducting state has a large number of degrees of freedom at low 
temperatures. So many, in fact, that it has completely swamped out the z scaling behavior. 
If we proceed to extremely low temperatures, five or six orders of magnitude below T c , we 
will eventually see the expected scaling. The low T c combined with the independence of z 
from T c has forced the effects of z to be pushed down to (presumably physically irrelevant) 
very low temperatures. 

There is a further reason why the superconducting state has so many degrees of freedom 
at low temperatures in these plots. This is the ground state entropy density of the normal 
state. According to the entropy balance, these degrees of freedom must reappear somewhere 
in the superconducting state. From the entropy balance formula (3) we see that the integral 
of cs(T)/T must equal the integral of cn(T)/T + stv(0)/T c . A nonzero sjv(O) and a low T c 
therefore push up the curve against which the superconducting phase degrees of freedom 
must balance. This curve is shown as a dashed yellow line in figure 2. As we recalled 
above, immediately below T c the superconducting phase must have a higher specific heat 
than the normal phase. However, in the presence of a ground state entropy density, this is 
not sufficient for the specific heat to need to decrease at low temperatures. In figure 2 we 
see that the specific heat of the superconducting phase is below the dashed yellow line at 
T c and is therefore forced to increase as the temperature is lowered, in order to make the 
entropy balance work. The ground state entropy of the normal state is therefore directly 
responsible for the temperature dependence of the specific heat in the superconducting 
state being markedly different to the prototypical BCS curve of figure 1. An increase in c/T 
as the temperature was lowered below T c in the ordered phase was one of the anomalous 
experimental features found in [2] , although no evidence was found there for a ground state 
entropy density. We will return to experiments in the discussion. 

Note that cjy/T is a constant in figure 2 not because the normal state is a Fermi liquid, 
it is not, but because the critical temperature is so low in this case that the specific heat of 
the normal state does not vary over the range < T < T c . 

Let us now increase the critical temperature by increasing 7. Figure 3 shows c/T against 
T / ji. We again take z = 1.5, 2, 3 but now with 7 = 2, 3, 5. The main feature of these plots 
is that as T c is increased, the anticipated low temperature scaling governed by z becomes 
visible. We have checked that the numerics reproduce (19). This is very much consistent 
with our discussion in the previous paragraphs. As T c is increased the system (i) has a 
wider range of temperature over which to balance the entropy and therefore the z scaling 
does not get swamped out and (ii) the effect of the normal state zero temperature entropy 
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Figure 3: Specific heat over temperature as a function of temperature for the normal (red) 
and superconducting (blue) phases below T c . This is for a model with a ground state entropy 
density sn(0). Three values of 7 and z are shown. The shaded region between the blue 
curve and the dashed yellow line, indicating cn/T + sjv(0)/T c , integrates to zero area. 
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density on the entropy balance is weakened according to sjv(0)/T c . 

For z < 2, the specific heat over temperature must still raise in the superconducting 
state as the temperature is lowered below T c , in order to cross the dashed line and achieve 
entropy balance, but then decreases at low temperatures as dictated by z. For z > 2, cs/T 
is monotonic. 

In this subsection we have seen how important features of the thermodynamics of the 
superconducting phase in these models can be understood by reference to the thermody- 
namics of the normal state. We saw that the nonzero sat(O) has a substantial impact on 
the specific heat in the superconducting phase, in some regimes swamping out the quantum 
critical physics. In the following subsection we turn to a quantum critical holographic model 
without a zero temperature entropy density in the normal state. 

4.3 Models with critical normal states with zn finite 

There are several ways in which to obtain a normal state with finite zn criticality. We use 
zn to denote the dynamical critical exponent of the normal state, in order to differentiate 
it from the z of the superconducting state that we discussed at some length in the previous 
subsection. The simplest setup seems to be, following [24, 25], to add a dilaton field to our 
previous action (27) 



For ease of comparison we will take the same potential and mass for the charged scalar 
as we did in (15) above. For simplicity we have chosen not to include a potential for the 
dilaton. Such a potential leads to interesting IR geometries [26], whose interplay with 
superconducting instabilities deserves study at some future point. Taking the same ansatz 
for the fields as in (20), with in addition the dilaton 3>(r), there are now five equations of 
motion given by 



C = 



1 

2/? 



R + 



6 

L2 



2 (V$) 2 j - -L^F^F"" - 1 (|V0 - iA4>\ 2 + mM 2 + V(\<f>\)) ■ 

(27) 




(28) 



(29) 



(30) 



(31) 



(32) 
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In the action (27) we have chosen the dilaton to be massless, and it therefore corresponds 
to a marginal operator in the dual field theory with dimension A = 3. We do not wish to 
turn on a coupling for this operator, and will therefore look for solutions to the equations 
of motion such that <I> ~ r 3 near the asymptotically AdS boundary at r — > 0. With this 
extra tuning to perform, the numerics otherwise proceed as before. 

In these theories the normal state, with (f> = 0, also has to be found numerically. How- 
ever, the low temperature, near horizon normal state geometries may be found analytically 
[24, 25]. In particular, the near horizon zero temperature metric has a Lifshitz form, and 
consequently one finds a specific heat scaling of 

c N ~ T 2 I ZN , z N = 1 + ^ • (33) 

a 2 

More specifically, the near horizon, r — > oo, normal state solution at zero temperature takes 
the form 

1 ,2 dt 2 dr 2 dx 2 + dy 2 . dt x 2, 

-j2 ds =- r 2 ZN + 9 L ^ r 2 + -2 ' A = ^^z N ^ $= a logr. (34) 

Where now = (1 + Zn)(2 + zn)/Q and h 2 L = (zn — 1)/{zn + 2). 

Finite zn means that these theories have many fewer degrees of freedom in the normal 
state at low temperature compared to the models of the previous section. According to the 
entropy balance we might, therefore, anticipate that T c will be lower. This expectation is 
too naive in general, as the entropy balance also depends on the low temperature degrees 
of freedom in the superconducting state, which we have not computed yet. However, it 
turns out to be true here. Figure 4 below shows T c //j> as a function of 7 for zn = 1.5, 2 
and 3, as well as for the locally critical model of the previous section. As previously, the 
critical temperature is independent of the quartic coupling u, as it only depends on the 
linearized equations of motion for the charged scalar. The plot is generated by looking for 
the temperature at which normalizable solutions to the linearized equation for the charged 
scalar field about the dilaton charged black hole background first exist. In the plot we see 
firstly that at fixed 7, larger z^ gives a larger T c . We will see shortly that this is loosely 
analogous to larger gap giving larger T c in the BCS case. Secondly, we see that as 7 becomes 
large, T c /fi tends to a constant value that is independent of zn- This second observation 
can be understood mathematically from the fact that the large 7 limit is the probe limit, 
in which the backreaction of the matter fields on the metric becomes negligible [15, 27]. 
In the probe limit, T c is determined uniquely by the Abelian-Higgs theory in a neutral 
Schwarzchild-AdS background to be T c j ' jj, 0.058 [15], in agreement with figure 4. 
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Figure 4: T c j \x as a function of 7 for, from left to right, z^ = 00, 3, 2 and 1.5. 

It is interesting that, in this class of models, T c j ' fi is bounded above by a rather small 
value. In the locally critical model, the critical temperature was obtained numerically as 
a function of 7 and the mass of the charged scalar field in [17], but the quantity plotted 
there was jT c /fx, so the universal limiting behavior (at fixed mass) was less immediate. 
While there is some ambiguity in the normalization of the charge density, it is natural to 
use (as emphasized in [17]) units in which charges are integrally quantized. This is the 
normalization we are using here. 

In locally critical models, the existence of a superconducting instability at sufficiently 
low temperatures can be diagnosed from the behavior of the charged scalar field in the 
near horizon extremal AdS2 X M 2 spacetime [16, 17, 28]. The charge of the scalar shifts its 
effective mass squared by the negative amount g tt AtAt [14], which can help push the effective 
mass squared below the ACLS2 Breitenlohner-Freedman bound, leading to a superconducting 
instability. In fact, the mass we have taken, m 2 L 2 = —2, already gives an instability of 
AdS2 at 7 = where the scalar is neutral [16]. We can see this in the zn = 00 curve 
in figure 4. In the Einstein-Maxwell-dilaton normal states being considered here, the IR 
effect of the charge of the scalar field is weaker. It is clear from the extremal near horizon 
solution (34) that in these models, at zero temperature, g tt AfAt — > in the near horizon 
limit r — > 00, and therefore the fact that the field is charged does not influence its IR mass 
squared. For the choice of mass m 2 L 2 = —2 we can look for an IR instability by solving the 
linearized equation of motion for the scalar in the near horizon Lifshitz background (34). 
The result is (independently of 7) 

^ r l+z N /2±^I%/(2V3) ( 35 ) 
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A complex exponent is associated to an instability [17, 22, 21, 9]. Therefore we see that 
there is an IR instability of the normal state against condensation of the scalar 4> whenever 
zn > 2. This is again consistent with the naive notion that larger zn favors instability. For 
zn < 2 the instabilities shown in figure (34) are not directly associated to the IR quantum 
critical regime, but depend on the full geometry. It is perhaps worthy of note that zjy = 2 
has the same number of low energy degrees of freedom as a Fermi liquid, as measured by 
the specific heat, while zn < 2 has fewer. 

Before moving on to numerical results for these models, we can determine analytically 
the low temperature behavior of the superconducting phase. As in the previous model, the 
low temperature behavior will be controlled by the near horizon geometry of the extremal, 
zero temperature, solution. Previously we found a Lifshitz geometry (16). The case with the 
dilaton is significantly more complicated because the dilaton and the charged scalar couple 
very differently to the Maxwell field. In particular, a purely Lifshitz type behavior combining 
the logarithmic running of the dilaton in (34) and the constant value of the charged scalar 
in (16) is not possible because the behavior of the Maxwell field is incompatible. Instead 
we take our cue from an insight in [21], where a similar technical challenge was faced for 
charged fermions. It is possible to find solutions where in the far IR the charged sector, i.e. 
the Maxwell field and charged scalar, are treated as a perturbation of a Lorentz-invariant 
IR spacetime. Thus we will find an emergent Lorentz invariant scaling with z = 1 in the 
superconducting state at low temperatures. 

One can solve the equations of motion (28) with the following expansion as r — > oo: 

/ = ^( 1 +/i r3 "° + '") • W 

' ~ 5T5=?( 1 + »'~ + ->. (37) 

h = h ir (1 - a ^ 2 + • • • , (38) 
<t> = \[\ + Pir^ a + --- , (39) 



2a S 48 7 2 



$ = -^iog v - ;r 2 ; +^ 3 - a +---- (40) 



The exponent a is not fixed; in order for this expansion to be valid we must have a > 3. 
The constants {/i, g±, h±,pi, P±} are easily determined but we will not need them. The 
above expansions describe an emergent IR AdS& spacetime in the superconducting state 
at zero temperature. This a familiar phenomenon, e.g. [29, 30, 31]. Our experience here 
suggests that the emergence of an AdS± may be generic in models with both dilatons and 
charged scalar fields. The incompatible couplings of the two scalars to the Maxwell field 
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effectively force the Maxwell field to go to zero in the IR. This indicates that the charge 
density operator is becoming irrelevant, and hence Lorentz symmetry re-emerges. 

The emergent Lorentz symmetry implies that the specific heat at low temperature in 
the superconducting state will go like 

c s ~T 2 . (41) 

These models therefore have in common with the previous locally critical normal state 
models that the critical modes are not entirely gapped in the superconducting state. An 
interesting difference is the following: In the model with a dilaton, the normal state has a 
tunable zn while the superconducting state has z = 1 fixed. The locally critical model had 
zn = oo fixed but z tunable. More interesting, however, is that in both cases we have 

z < z N . (42) 

While, like most of our observations, this inequality does not follow as a necessity from 
entropy balance, it is nicely in agreement with a simple-minded application of the need to 
balance the degrees of freedom; the superconducting phase has more degrees of freedom than 
the normal phase at higher temperatures, just below T c , and so should have fewer degrees 
of freedom at lower temperatures. A smaller z indeed implies fewer low energy degrees 
of freedom according to (19). It would be interesting to see how robustly the inequality 
(42) holds in holographic models in which both the normal and superconducting dynamical 
critical exponents can be tuned. Models in which the normal state Lifshitz geometry is 
achieved via charged fermions rather than dilatonic scalars [32, 33] may provide a framework 
to address this question. 

In figure 5 we plot the specific heat over temperature for the normal and superconducting 
states for temperatures below the critical temperature. We chose the three values of zn = 
1.5,2,3 in order to illustrate cases in which the specific heat over temperature goes to 
zero, a constant and infinity at low temperatures, according to (33), respectively. These 
behaviors are clearly visible in the plots. We have already noted that at fixed 7, the critical 
temperature increased with zn- In figure 5 we have chosen 7 so that the three values of 
T c /fx shown are the same for each value of zjy and furthermore comparable to the critical 
temperatures we plotted for the locally critical case. 

The anticipated linear behavior of C5/T at low temperatures is also clearly visible in 
the plots. We have checked that the numerics reproduce in detail the low temperature 
behavior of equations (36) - (40). In fact, the curves for cs/T are essentially featureless 
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Figure 5: Specific heat over temperature as a function of temperature for the normal (red) 
and superconducting (blue) phases below T c . These plots are for models without a ground 
state entropy density. Three values of T c //x and zm are shown. The parameter 7 has been 
chosen in each case to make the three values of T C J ' \i equal. The shaded region between the 
blue and red curves integrates to zero area. 
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and decrease monotonically from their values at T c . These plots are closer to the BCS 
curve of figure 1 than the locally critical models. The main difference with the BCS case, 
beyond the fact that the ordered state is not gapped here, is that the shape of the curves 
changes as T c is increased; the jump in specific heat at the phase transition decreases with 
increasing T c . More generally, the two curves come closer as T c is increased, consistent with 
the fact, mentioned above, that in the large 7 limit the charged sector of theory becomes 
an increasingly negligible 'probe' sector that does not backreact on the neutral degrees of 
freedom. One way to phrase this is that 7 controls the static charge susceptibility of the 
theory, which is becoming small. 

5 Discussion 

In this paper we have used the notion of entropy balance to relate qualitative features of the 
normal and superconducting phases in various models of holographic superconductivity. As 
part of this investigation we have introduced new holographic superconductors in which the 
critical normal state from which the superconductivity emerges has finite dynamical critical 
exponent zn- While entropy balance itself is a rigorous statement, we have emphasized 
that the nature of the relations we ascribe to it is not one of logically necessary connections, 
but rather of tendential relations that one might 'naively' but reasonably expect to hold 
more often than not. The main observations connected to entropy balance that we have 
made for the models considered in this paper are that (a) z < zn, (b) larger zn leads 
to a larger T c if other parameters are kept fixed, (c) cs/T starts decreasing below T c in 
models without a ground state entropy, but increases in models with a ground state entropy 
density, (d) in models with a ground state entropy density and a low critical temperature, 
critical scaling in the superconducting state is pushed down to very low temperatures. The 
simplest way in which these statements can fail to hold is if the jump in specific heat at the 
critical temperature varies in such a way that it cancels out the effects of low temperature 
criticality in the entropy balance formula. It is of interest perhaps to understand more 
deeply the physics controlling this jump. 

Our observations cannot be applied directly to either the nematic or ferromagnetic sys- 
tems that we cited in the introduction as motivation for this work. The ferromagnetic 
systems have first order phase transitions. The nematic phase in Sr3Ru207 does have a 
second order transition, but the ordered phase develops out of an anomalous Fermi liquid 
normal state in which cn/T is logarithmically increasing as the temperature is lowered. If 
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this logarithm is extrapolated down to zero temperature, then it is able to precisely en- 
tropically balance the anomalous increase in the specific heat of the ordered phase, without 
recourse to a ground state entropy density in the normal phase [2]. This is ultimately tied 
to the fact that nematicity following from a Pomeranchuk instability, unlike superconduc- 
tivity, leaves the existence of a Fermi surface intact and therefore allows many more low 
temperature degrees of freedom in the ordered phase. An increasing specific heat below T c 
in a superconductor would be more dramatic. 

Heavy fermion superconductors are another class of systems with unconventional nor- 
mal and superconducting states, in which entropy balance has been explored experimentally 
[36, 37]. The heavy fermion materials exhibit power law specific heat at low temperatures in 
the superconducting state, and thus, like the models we have considered, have more degrees 
of freedom at low temperature than gapped BCS superconductors. Unlike the holographic 
models considered in this paper, which have scalar condensates, the extra degrees of free- 
dom in the heavy fermions are naturally associated with nodal excitations of the non-s-wave 
Cooper pairing. Somewhat similarly to S^R^Oy, these extra degrees of freedom are en- 
tropically balanced out with additional degrees of freedom in the non- Fermi liquid normal 
state [36, 37]. 3 In CeCoIn 5 there appears, in addition to the power law in temperature, 
a rather spectacular surplus of very low temperature degrees of freedom in the supercon- 
ducting phase [38]. The entropic connection between unconventional normal states and 
unconventional superconductivity in these materials is precisely the type of physics we have 
attempted to get a computational handle on in this paper. It may be possible to make a 
tighter connection to the materials if the holographic emergence of superconductivity from 
a fractionalized phase can be related to, for instance, spin density wave physics, in the spirit 
of e.g. [11]. 

Specific heat measurements also exist for the cuprates and show some similar features, 
e.g. [39]. In this case, extracting sharp information from the measurements is complicated 
by the need to subtract a phonon background, due to the higher temperatures involved. 

Beyond the models we have considered here, there are other holographic contexts where 
the interplay of superconductivity and quantum criticality may be usefully studied. Quan- 
tum critical normal states with zn finite can be constructed using 'electron stars' [32, 33] 
rather than dilatonic fields. These are qualitatively different in that the asymptotic electric 

3 The way in which entropy balance works out in these materials suggests that T c would be increased if 
the non-s-wave condensates were to emerge from a strict Fermi liquid normal state, without a logarithmic 
enhancement of the low temperature degrees of freedom. 
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flux is sourced by charged fermions in the bulk rather than emanating from behind a hori- 
zon [9]. Thus any superconducting instability in this case is not coincident with exit from 
a fractionalized phase, but corresponds to the re- arrangement of confined charge- carrying 
'quark' degrees of freedom from fermionic to bosonic 'mesons' [9]. This phenomenon has 
been studied for neutral scalar instabilities in [22], but not yet for charged scalars. 

A further way to probe the interrelation between deconfined critical degrees of freedom 
and superconductivity is to incorporate the physics of confinement. Holographic models 
have been developed that exhibit confining as well as superconducting phases [34, 35]. It 
is of interest to extend our computations to include these effects. In particular, following a 
possibly overly naive application of entropy balance, one would like to confirm the extent to 
which gapping out the critical modes in the normal state causes a reduction in the critical 
temperature. 
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